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ABSTRACT 



Context. The determination of horizontal velocity fields at the solar surface is crucial to understanding the dynamics and 
magnetism of the convection zone of the sun. These measurements can be done by tracking granules. 

Aims. Tracking granules from ground-based observations, however, suffers from the Earth's atmospheric turbulence, which 
induces image distortion. The focus of this paper is to evaluate the influence of this noise on the maps of velocity fields. 
Methods. We use the coherent structure tracking algorithm developed recently and apply it to two independent series of images 
that contain the same solar signal. 

Results. We first show that a k — u filtering of the times series of images is highly recommended as a pre-processing to decrease 
the noise, while, in contrast, using destretching should be avoided. We also demonstrate that the lifetime of granules has a strong 
influence on the error bars of velocities and that a threshold on the lifetime should be imposed to minimize errors. Finally, 
although solar flow patterns are easily recognizable and image quality is very good, it turns out that a time sampling of two 
images every 21 s is not frequent enough, since image distortion still pollutes velocity fields at a 30% level on the 2500 km scale, 
i.e. the scale on which granules start to behave like passive scalars. 

Conclusions. The coherent structure tracking algorithm is a useful tool for noise control on the measurement of surface horizontal 
solar velocity fields when at least two independent series are available. 
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1. Introduction 

Movies of the solar surface show that it is a place of intense 
turbulent fluid flows where three major scales (granula- 
tion, mesogranulation, and supergranulation) have been 
pointed out. In order to better understand the underlying 
dynamics, it is crucial to be able to measure the velocity 
fields. As far as horizontal flows are concerned, the basic 
techniques have relied on m easuring the di s placem ent of 
granules that, as shown by iRieutord et ah ( 2001 ). trace 
the fluid flows on scales larger than 2.5 Mm. 

Two algorithms have been devised to transform a 
time sequence of images into a sequence of horizontal 
velocity fields maps. These ar e the LCT algorithm (i.e . 
local correlation tracking, see iNovember &: Simon 1988) 
and the CST algorit hm (coherent structure tracking, see 
Rieutord et al.ll2007l hereafter referred to as Paper I). 
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When the surface velocity field is known, one is usually 
interested in identifying/following the dynamical struc- 
tures of the flow - like vortices, upwellings, or down- 
wellings. The identification of these structures demands, 
however, computing the velocity gradients like the diver- 
gence or the vorticity. In Paper I, it has been pointed 
out that such quantities are very sensitive to the noise 
induced by terrestrial atmospheric distortion, since they 
are derivatives of the velocity field. The use of the veloc- 
ity field to describe the dynamics of the solar surface thus 
needs to be complemented by an error analysis that both 
evaluates the significance of the observed dynamical fea- 
tures and gives a way to eliminate, or at least reduce, the 
impact of errors. 

The aim of this paper is to analyse the conse- 
quences of errors in the final result of velocity, vortic- 
ity, and divergence fields. As mentioned above, the main 
source of errors comes from the distortion of images in- 
duced by the Earth's atmospheric turbulence. The case 
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Fig. 1. View of the region used for the tests. X and Y 
scales are in arcsec. 



of errors or, equivalently, the precision of measurements 
has already been d i scuss ed for the LCT algorithm by 



November fc Simonl 



who mentioned errors the or- 



der of 20 m/s on the velocity field. Further work by 
Simon et al. I (|l995l) showed that this precisi on was cer- 
tainly largely overestimated. More recently, iPotts et"al 



( 20031 ) investigated the case of interpolation errors, as- 
sociated with the LCT algorithm, which also spoil the 
final result. Here, we focus on the CST algorithm and try 
to give a neat picture of the influence of the Earth's at- 
mospheric distortion on the measurement of the velocity 
fields on different scales. 

We organised the paper as follows. Using two indepen- 
dent series of images of the solar surface, we first evaluate 
the noise induced on the positions of the granules and how 
image pre-processing can reduce it. We then focus on the 
way the noise influences the final velocity fields on differ- 
ent scales and analyse its propagation up to the curl and 
divergence maps. Conclusions and outlooks follow. 

2. Observations 

We use a time series of images obtained on June 5, 
1993 at the SVST (Swedish Vacuum Solar Telescope), 
Observatorio Roque de los Muchachos, La Palma 
(data courtesy of P. Brandt, see also Simon et al.l 
19941) . This da ta set has al r eady bee n studied by 



many a uthors dSobotka et al. 1997allb ; Sobotka et al 



1999allbl: iDorotovic et all 120021: iGetling fe Brandt] l2002t 
Rou dier et al .1 120031 ). The original series consists of 1868 
image pairs of size 1310x970 pixels taken at A = 468 ± 
5 nm. The time between two pairs of images is close to 21 
seconds. Images from a pair are separated by a few sec- 
onds (3 s on average and always less than 14 s). The pixel 
size is 0'.'125, and the spatial resolution is near the diffrac- 
tion limit 0"25. The field of view is 2f7x2.'0. However, the 
instrument leads to a rotation of the field of view, and 
the area observed on the Sun at the beginning of the time 
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Fig. 2. Indices of image quality. Contrast (solid line) 
and distortion amplitude (dashed line, scale on right). 
Distortion is the rms displacement, evaluated with local 
correlation tracking, between two images of a pair. 
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Fig. 3. The normalized intensity contrast versus the rms 
displacement of granules between the two images of each 
pair. 



scries is different from the one at the end. The rotation 
centre is located on a pore at pixel coordinates (590, 102) 
(not in the field-of-view used here). 

For our investigations described below, we used a sub- 
sample of 210 pairs of frames (images 841 to 1262), cover- 
ing ~ 77 minutes. In Sect. 6, slightly more data have been 
used (images 800 to 1298 corresponding to a ~ 87 min 
sequence). As shown in Fig. [TJ we extract a region of 
401 x 401 pixels centred on a (magnetically) quiet zone. 
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Fig. 4. Mean distance between granule centres of grav- 
ity for each image pair versus the mean displacement due 
to atmospheric turbulence estimated by local correlation 
(one point per image pair). Both displacements are in pix- 
els. The solid line shows the line of equation y — x and 
the dotted line is a linear fit over the points. 



The main advantage of this data set is that it con- 
tains two independent sequences of images that can be 
considered as representing the same solar signal. The only 
difference comes from the Earth's atmospheric distortion, 
the effects of which can thus be analysed. Moreover, as 
shown by Fig. [21 image quality varies during the sequence. 
Although the contrast remains almost constant, we see 
that the amplitude of distortion increases after t=45 min. 
We thus have at our disposal a "good" sequence where 
mean distortion is about 0.5 pixel and a "bad" sequence, 
where mean distortion reaches an amplitude of 2 pixels. In 
the following when we refer to "good" data, we mean the 
first 45 min, while "bad" data will designate the remaining 
sequence. 

3. Sources of errors in CST 

In the CST method, errors are introduced through the 
segmentation (i.e. through determination of the granule 
positions) and through interpolation done in order to ob- 
tain a velocity field sampled over a regularly-spaced grid. 
The precision at which the velocity field can be measured 
thus depends on: 

— the precision at which a granule position is determined, 

— the duration of granule tracking, 

— the temporal resolution (time interval between velocity 
maps), 

— the spatial resolution of the grid used to sample the 
velocity field. 

This list shows that the CST algorithm allows us to iden- 
tify all the crucial steps and thus to follow the propagation 
of errors from the beginning to the end of the computa- 
tion. 



4. Standard deviation of granule positions 

The first step in the error analysis is to estimate the er- 
ror on the granule position^. We first apply it to the raw 
data, and then use it to estimate the performances of pre- 
processing steps that can be performed before granule seg- 
mentation, namely on the image sequences. 



4.1. Method 

To estimate the standard deviation on granule positions, 
we use pairs of images. Because of their quasi-simultaneity, 
they represent the same solar surface, and all the differ- 
ences between them come from atmospheric distortions, 
instrumental effects, or processes applied on the data. 

After the segmentation step, the algorithm identifies 
granules present in both images of the pair. We then sim- 
ply measure the distance between two identical granules. 
More precisely, we consider that the position of a granule is 
controlled by two random variables (xi,yi) for which we 
have two realisations. We thus construct two other ran- 
dom variables (5x l: Syi) = (x,,i - x h2 ,yi,i - 2/4,2), whose 
variance is just twice that of (xi, yi). Thus with a pair of 
images we can estimate the mean error over the field of 
view in each direction x and y by 
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and the mean displacement Sr = 

These quantities are clearly estimates of the im- 
age quality through atmospheric distortion. Interestingly 
enough, we compared this estimate with the intensity con- 
trast usually used to indicate image quality. As shown in 
Fig- El the correlation between the two indicators is rather 
poor, meaning that they are largely decoupled although 
they both come from atmospheric turbulence! We inter- 
pret this result, tentatively, as showing that different lay- 
ers of the Earth's atmosphere control the contrast and the 
distortion. 

Furthermore, we compared (see Fig. [4]) the displace- 
ment of granules with the displacement field derived from 
a local correlation tracking (using an FWHM of 20 pix- 
els) between the two images of a pair. As expected, the 
correlation is much better (~ 0.96). 

4.2. Application to pre-processing 

One of the interesting applications of the error estimate 
on the granule position is to allow the evaluation of the 
effectiveness of various pre-processings applied to the im- 
ages before estimating the granule position. A measure of 
<j x before and after a given pre-processing indicates its 
performance. We have applied this approach to the k — lu 
filtering and to the destretching. 



By granule position we mean the position of the centre of 
gravity of the granule in the segmented image. 
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Fig. 5. Estimation of the error a x on the determination 
of the centre of gravity position in the x direction as a 
function of time with no k — lu filtering (solid line) and 
with the filtering (dotted line). 







Good data 








a x Oy 


N g 


No 


processing 


0.696 ±0.081 0.701 ±0.084 


1077 ± 91 




k — LU 


0.416 ±0.037 0.418 ±0.028 


1096 ± 34 






Bad data 








o x Oy 


N g 


No 


processing 


0.694 ±0.075 0.727 ±0.079 


803 ± 150 




k — LU 


0.470 ±0.039 0.527 ±0.045 


972 ± 89 



Table 1. Comparison of the performances obtained on 
time series with high and low quality images with and 
without k — lu filtering. 



4.2.1. k — lu filtering 



The k — lu filtering (see for example lTitle et al.l . ll987T ) acts 
both spatially and temporally. In practice, it is a thresh- 
olding in the Fourier space. The contributions of all fre- 
quencies corresponding to a phase velocity higher than the 
threshold are eliminated. 

Let us note that the data we use are irregularly sam- 
pled. The time step between two pairs of images is about 
21 seconds, within a few seconds (the maximum deviation 
with respect to a periodic sampling is 15 s, while the av- 
er age deviation is 6 . 3 s). T he same data set has been used 
by iDorotovic et ah ( 20021 ). who interpolated them to get 
a regularly spaced time series. We consider that this is an 
unnecessary refinement since, within a few seconds, the 
solar signal does not change. Hence, with a regular time 
step, images would differ from ours by just a different re- 
alisation of the Earth's atmospheric turbulence. 

Figure[5]shows the average pixel error on the x-position 
of the granules for data pre-processed with a fc — lu filtering 
and data without pre-processing. We used a threshold for 
the phase velocity of 4 km/s. The figure shows that the 
performance is quite improved over the whole time series 
by the use of this filtering. The error for the y-component 
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Fig. 6. An example of a map of displacements due to the 
Earth's atmospheric distortion as sampled by granules. 
Regions without data have too low an image quality for a 
granule to be identified. X and Y are in arcsec. 



is within an order of magnitude of the error for the x- 
component. 

One may wonder why the increasing distortion seen 
after t=45 min in Fig. [2] does not appear in Fig. \E\ This 
comes from the way the error is measured. Indeed, only 
granules close enough in a given pair of images are kept. 
The distortion map in Fig. [5] shows the lack of measure- 
ments in some regions: image quality is not good enough 
for granules to be identified from one frame to the other. 
With this method it is clear that high-amplitude distor- 
tion does not show up with increased error in granule po- 
sition, but instead appears in a reduction of the number 
of "valid granules" . In Fig. [7] we see the loss of granules 
when conditions deteriorate and how k — lu filtering im- 
proves the situation. Numbers in Table Q] also illustrate 
this process. The errors only slightly increase from the 
"good" to the "bad" data, but the number of granules, 
and in other words the number of measuring points, is re- 
duced by 20% without filtering, while the loss is 10% with 
filtering. 

4.2.2. Destretching 

This method, introduced by I November &: Simon (|l988lh 
is based on the same principles as the LCT. It uses a local 
correlation scheme to determine the displacements with 
respect to a reference image. From these displacements, 
images are stretched by interpolation in order to be the 
closest possible to the reference image. This method aims 
at correcting the effects of atmospheric turbulence and 
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Fig. 9. Efficiency of destretching: the dotted line shows 
the difference between two simultaneous images when one 
has been destretched with the other as reference. 



Fig. 7. Comparison of the number of granules for each 
image for raw data (dotted line) and for k — u> filtered 
data (solid line). The number of granules is defined as 
the number common to both images taken 'simultane- 
ously'. We note that filtering attenuates the fluctuations 
and when atmospheric distortion increases, filtering re- 
duces the losses of granules. 



1.2 
1.0 



=3 0.8 



— 0.6 

£ 0.4 

0.2 
0.0 







20 40 

time (min) 



60 



Fig. 8. Same as Fig. [5] but for destretching; the solid line 
is for unprocessed data, while dotted lines show the error 
{(Jx) with the destretched sequence. 



allows us to compare images that are then plagued with 
the same distortion. 

The pre-processing we used applies 4 successive correc- 
tions based on local correlation tracking with the following 
parameters: an FWHM of 31 pixels with a 62 pixel step, 
an FWHM of 62 pixels and a 31-pixel step, a 32 FWHM 
and a 16-pixel step, and then an FWHM of 20 pixels and a 
10-pixel step. This combination has been established em- 



pirically and seems to give the best results. We remind 
that 62 pixels correspond to 7.75" and 10 pixels to 1.25". 

Figure [8] compares the error estimate on the granule 
positions in the x-direction between the destretched and 
raw time series. Note that when image quality is good 
(t < 45 min) errors remain of the same order of mag- 
nitude as in raw data but fluctuations are less impor- 
tant, while when image quality decreases (t > 45 min), 
errors increase. Quite clearly the destretched series does 
not compare favourably with the raw one. On the con- 
trary, destretching seems to worsen the results when the 
image quality is slightly degraded. 

Actually, the poor performance of destretching on er- 
ror reduction could be anticipated from the result dis- 
played in Fig. There we plot the mean distance be- 
tween the granules in two simultaneous plates when one 
plate has been destretched to the other. If destretching 
were perfect, the error would vanish. Clearly, this is not 
the case: in the sequence with 'intense' atmospheric tur- 
bulence the error is not reduced at all, while during the 
good sequence, a small factor 1.5 is gained. Therefore we 
interpret the error increase generated by destretching as 
evidence that the destretching process decorrelates from 
the true displacement of granules and thus introduces a 
new random variable whose dispersion adds to the origi- 
nal signal as shown by the factor taken by the error. 
This decorrelation may probably come from a change in 
scale of the distortion motions that is no longer matched 
by the destretching process optimized for the first frames 
of the sequence. If this interpretation is correct, the de- 
stretching process would need a readjustment of the local 
correlation tracking parameters from time to time, mak- 
ing the whole processing very costly computationally. The 
only good point introduced by destretching is the reduc- 
tion of the fluctuations. 
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Table 2. Rms fluctuations of the position of granules on a 
typical pair of images shown for different resolutions with 
the multi-resolution analysis and with and without k — ui 
filtering, as are given in image pixel. The 1'.'25 resolution 
corresponds to unfiltered data. 



4.2.3. Distortion noise on different scales 

We have shown in Paper I that multi-resolution analy- 
sis was an interesting tool for the determination of flow 
structures. We may thus wonder how the distortion noise 
affects the different scales and therefore resists to wavelet 
filtering. Using maps such as the one in Fig. [6] we compute 
the different components of a multi-resolution decompo- 
sition as will be used in the analysis of the velocity field. 
Quantitatively we show in Table[2]the values of the distor- 
tion amplitude on different scales. This table again shows 
the importance of fc — uj filtering in the reduction of errors, 
especially on large scales. 

4.2.4. Conclusion 

To conclude this section, two points should be underlined. 
On the one hand, destretching is certainly an unnecessary 
complication whose only positive effect is to reduce the 
noise fluctuations when correctly tuned to the distortion 
scales; otherwise, it is likely to double the variance of the 
signal. On the other hand, k — uj filtering appears as the 
required pre-processing whose effect on noise reduction is 
clear. We surmise that, when using time-sequences with a 
higher time-sampling, the noise reduction by fc— lu filtering 
will be even more efficient. 



5. Propagation of errors from granule positions to 
the velocities 

After examining the estimation of errors on the granule 
position in the previous section, we now study the propa- 
gation of these errors to the velocity of the granules, after 
averaging in time and over a given spatial range. 
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Fig. 10. Error on the x-component of the velocity as a 
function of the threshold on the granule lifetime. The 
different lines are for different resolutions of the multi- 
resolution decomposition; the solid line shows the raw case 
(no wavelet filtering). The optimal threshold clearly ap- 
pears around 3 min. 
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Fig. 11. Average distance between the granule centre of 
gravity for a pair of images as a function of the granule 
size. 



5.1. Error propagation 

As seen above, the first step of the velocity computation 
in the CST algorithm is to determine of the position of 
granules on each image. The duration of the tracking of 
the fc th granule is given by Atk = i n/ (fc) — ^n;(fc)' where 
rii(k) is the image where the granule k appears and nf(k) 
the image where it disappears. 

Over an image, granules are much less dense than pix- 
els, making the velocity field sampled on a much coarser 



grid whose elements have size 5. The velocity at a grid co- 
ordinate (x, y) is assumed to be the average of the velocity 
of granules whose average coordinates belong to the do- 
main D around (x, y) (i.e. in [x— 5/2, x+S/2], [y— 5/2, y + 
5/2] ). The spatial resolution of the velocity field is given 
by the mesh size 5. The x-component of the average ve- 
locity in D is given by 
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where AT is the number of trajectories falling in D. The 
uncertainty on the value of V x is then 



1 

k£D 



k,n f 



If we assume that the dispersion on the centre of grav- 
ity remains the same for the whole time series (<J X ) and 
that the time interval Atfc is the same for all granules (i.e. 
all granules have the same lifetime At), then the expres- 
sion of the dispersion of the average velocity in x is 



1 



2a? 



2a 2 



N 2 ^ At 2 

i£D 



NAt 2 ' 

In this case, the error on the velocity varies like 



SV = 



\/2<T x 



NAt 



(1) 



We see that precise velocity values need many granules 
in a grid element and a long time interval. In other words, 
errors are less if a coarse resolution in space and time 
is used. A trade off must be found (see below), but it 
is clear that this technique will be more appropriate to 
slowly evolving large-scale flows than to rapidly varying 
small-scale ones. 
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5.2. Influence of the granule lifetime 

However , all granules do not hav e the same lifetime (see for 
example iHirzberger et al.lll999l or Paper I). In this case, 
it may be more appropriate to select only a subsample of 
the granules. For example, in the simple case where half 
of the A~ granules have a lifetime Aat and the other half 
a lifetime Abt = 2Aa t , the error on the velocity is 



'a+b 




iNAaf 



If we take only the ^ granules with the largest lifetime 



Ab t into account, we obtain 



ieD,Ati=Ab t 
T 2 



4 Na 2 x 
N 2 4A a 2 
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NAa 2 



< a. 
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Therefore, the errors are smaller when considering only 
the granules with the longest lifetime. 
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Fig. 12. Two views of the divergence field at different 
resolution. The difference between the left and right flow 
fields is the noise introduced by the Earth's atmosphere. 
In the first row no filtering has been applied and common 
features are barely identifiable. The following rows show 
filtered data according to multi-resolution representation 
(see Paper I) with a resolution divided by 2 from one row 
to the next. The mesh size is 10 pixels (1'.'25) and the flow 
is an average over ~ 1.5hr. X and Y are in arcsec. 



This result shows that some selection of "valid gran- 
ules" may improve the quality of the velocity field. We 
thus impose a lower threshold on the granule lifetime to 
eliminate short-lived structures. This threshold has to be 
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Fig. 13. Same as in Fig. [TO] but for the z-component of 
the curl field. 

determined empirically and depends on the filling factor of 
the grid. (We need to avoid grid points with no data, see 
Paper I.) This strategy may be improved by computing 
the average velocity for various thresholds, starting with 
the highest one. Since granules with long lifetimes are rare 
the number N is small and the error is large. By reducing 
the threshold, the number of granules increases rapidly 
and the error decreases. For some optimal threshold, the 
error ceases to decrease as the increase in granule num- 
ber no longer compensates for the diminishing value of At 
(see Eq. [1} . This process is illustrated in Fig. [TO] where 
we show the dispersion of velocity differences between the 



two independent data sets as a function of the threshold 
on the granule lifetime. Clearly, for these data the optimal 
threshold is around 3 min. 

5.3. Influence of the granule size 

It is also interesting to study the influence of the granule 
size on the error for the granule centre of gravity. We would 
expect a smaller error in the case of large granules, as 
these are defined by a larger number of pixels. A selection 
on the granule size could then also improve the precision 
on the velocity field. Figure [IT] shows an estimation of 
the error on the granule position for various granule size 
intervals. The precision increases only slightly with the 
granule size, and the improvement is at most 0.05 pixels. 
This means that the error on the velocity will not change 
much either. This small improvement in the error is too 
weak to compensate for the increase in the error on the 
velocity due to the decrease in the number of granules as 
their size increases. 

6. Error propagation to the final maps 

Following the algorithm described in Paper I (but see also 
Sect. 5.1), we now compute the velocity fields of the two 
independent time series at a different resolution (see Paper 
I). A comparison between the two fields shows the influ- 
ence of the distortion induced by the Earth's atmosphere 
and the efficiency of MRA in revealing the flow patterns on 
the different scales. The results are illustrated by Figs. [TO] 
and !131 which show the horizontal divergences and the ver- 
tical component of vorticity, respectively. This flow field is 
the average velocity over a time lapse of At = 87 min in- 
cluding the preceding sequence and with no distinction be- 
tween the "bad" and "good" sequences. Actually, we first 
computed the velocity fields using the "good" sequence 
and then extended the computation to the whole set. We 
observed that the dispersion between the two independent 
sequences was still reduced when using the whole set of 
data showing that, for the determination of the large-scale 
mean flow, the observed degradation in the distortion is 
not influential. 

In Table [3] we quantitatively summarise the dispersion 
of the results as a function of the resolution for unfiltered 
images. This table has been obtained with the whole time 
series (87 min), using a velocity pixel of 1.25" and re- 
moving all granules with a lifetime less than 180 s. We 
computed the rms velocities, divergences and curls, of the 
maps issued from the multi-resolution decomposition. The 
"uncertainties" shown along the numbers give the am- 
plitude of the fluctuations generated by the Earth's at- 
mospheric distortion. This table shows the decreasing in- 
fluence of distortion with increasing scale. It also shows 
that, as expected, velocity gradients (divergence and the 
z-component of the vorticity) suffer much more from the 
noise and that the curl is certainly the quantity most sen- 
sitive to image quality. The values of the correlations be- 
tween the results issued from the two independent series 
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Resolution (pixel size) 


Vrms 


Cu 


(Div V) rms 


c d 


(Curl V) rms 


c c 


Turbulence 


1.25" 920 km 


700±695 


0.49 


7.1±7.0 


0.45 


5.7±7.0 


0.22 


880±224 


2.50" 1840 


478±342 


0.74 


3.8±2.8 


0.69 


2.9±3.0 


0.45 


893±122 


5.00" 3680 


340±172 


0.87 


1.6±0.9 


0.89 


1.2±0.9 


0.73 


890±94 


10.0" 7360 


205±92 


0.90 


0.58±0.24 


0.94 


0.41±0.24 


0.85 


896±65 



Table 3. Rms velocities, divergences and curls, of the maps issued from the multi-resolution analysis of the whole 
time series. The first row is raw data (no wavelet filtering) while the next rows shows the numbers issued from the 
wavelet-filtered maps. Velocities and turbulence are in m/s, divergence and curl in 10 _4 s _1 . 



Resolution (pixel size) 


Vrms 




(Div V) rms 


c d 


(Curl V) rms 


C c 


Turbulence 


1.25" 920 km 


597±330 


0.85 


6.8±4.2 


0.81 


5.0±4.2 


0.66 


680±193 


2.50" 1840 


466±162 


0.94 


4.0±1.7 


0.91 


2.7±1.7 


0.81 


692±105 


5.00" 3680 


350±80 


0.97 


1.8±0.4 


0.97 


1.2±0.4 


0.95 


689±78 


10.0" 7360 


209±40 


0.98 


0.63±0.1 


0.98 


0.42±0.1 


0.96 


694±53 



Table 4. Same as Table [3] but for k — ui filtered data. Note the reduction in the dispersion of the results. 



Resolution (pixel size) 


Vrms 


G v 


(Div V) rms 


Cd 


(Curl V) rms 


c c 


Turbulence 


07 875 644 km 


722±534 


0.73 


11±9.7 


0.64 


10±9.7 


0.54 


528±288 


1775 1288 


528±258 


0.88 


6.3±3.7 


0.83 


4.5±3.8 


0.65 


532±143 


3750 2576 


403±128 


0.95 


2.7±0.9 


0.94 


1.8±0.9 


0.87 


533±75 


7" 5152 


284±65 


0.97 


1.2±0.3 


0.97 


0.7±0.3 


0.94 


533±45 



Table 5. Same as Tabled] but for a smaller velocity pixel. 



(C„, Cd and C c ) quantitatively indicate the similarity of 
the fields. Here too, filtered fields are much better corre- 
lated, up to 90% on the velocity field. As shown below, 
this correlation is still improved when using k — lo filtered 
data. 

In this table, we also give a 'turbulent velocity', which 
is the mean rms dispersion of granule velocities. Indeed, 
in each velocity pixel we take the mean velocity of the 
granules falling in this very pixel; however, while comput- 
ing this mean velocity, we also have access to the disper- 
sion around this mean. This mean (over the whole field of 
view) dispersion represents the random motion of granules 
around their drift by large-scale flows. This table shows 
that, although this quantity suffers (also) from image qual- 
ity, it is almost independent of scale. This independence 
is expected since this quantity measures the proper mo- 
tion of granules and therefore their intrinsic kinetic energy, 
which should not vary from place to place. 

Table 2] gives the same quantities but for the k — us fil- 
tered sequence. We note the strong reduction of the noise, 
almost a factor 2, on all scales and the strong improvement 
of correlation between the results of the two independent 
series. Moreover, the dispersion of the velocities is reduced 
in the same proportion as displacement of granules (i.e. 
e95™^ ~ °'^ ). The comparison of the values of turbu- 
lence in Tab. [3] and 2] shows the influence of the Earth's 
atmospheric noise on the random motion of granules. 



Table [5] gives another view of this velocity field us- 
ing a smaller mesh size. Here we computed the velocity 
amplitudes that are traced by granules when they can be 
considered as passive scalars, i.e. on a scale larger than 
2500 km. Numbers show that on that scale the velocity 
field has an amplitude of 400 m/s and that such a mea- 
surement is still uncertain by 30%. This is about the same 
for the divergence, but it rises to 50% for the vorticity. The 
more intense fluctuations compared to Table 3] come from 
the smaller scales involved. The weaker "Turbulence" val- 
ues come from the empty bins which are more numerous. 

These results show that, when independent images se- 
ries are available, the use of the CST algorithm authorizes 
a tight control on the role of the noise induced by the 
Earth's atmospheric distortion. 

7. Conclusions 

In this paper we have tried to show the various influences 
of the noise generated by the Earth's atmosphere when 
one measures the horizontal velocity field at the surface of 
the sun through granule tracking. 

For this purpose we compared two time series of images 
of the solar surface separated by at most a few seconds; 
they thus represent the same solar signal but with a dif- 
ferent noise (essentially, distortion from the Earth's atmo- 
sphere) . We measured the position of the granules in the 
two time series and used these data to determine the am- 
plitude of the noise. We could thus test the pre-processing 
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that could first be applied to the images, namely a k—u> fil- 
tering and the destretching. The comparison between the 
pre-processed series and the raw one allowed us to evaluate 
the efficiency of the pre-processing. It turns out that k — u> 
filtering significantly reduces the noise while destretching, 
even if reducing the fluctuations, not only cannot reduce 
it, but amplifies it when distortion increases. 

We also found that the CST algorithm that tracks 
coherent structures (essentially granules) could easily go 
across sequences of images with decreased image quality. 
This is because all granules are not affected evenly by dis- 
tortion and granules whose trajectories are too perturbed 
by atmospheric noise are eliminated; they thus do not in- 
put noise in the final velocity field. This is clearly a feature 
that algorithms based on local correlations of images can- 
not authorize. 

We also studied how the distortion noise affecting 
granules positions introduces some noise into the inter- 
polated velocity fields. We thus showed that the lifetime 
of granules was an important parameter and that short- 
lived granules should be eliminated. The trade-off between 
granule number (the more numerous the granules, the bet- 
ter sampled the velocity field) and noise intensity seems to 
be, for our sequence, around 3 min. It is clear that, if the 
atmospheric noise is more intense, this threshold should 
be increased. With the decomposition of the velocity on 
the Daubechies wavelets (see Paper I), we could evaluate 
the impact of the noise on different scales and show that 
errors on velocities decrease with increasing scale, as ex- 
pected. More precisely, we could show that on a scale of 
2500 km, i.e. on the scale where granule m otions trace 
the large-scale flows ( Rieutord et all l200ll ). the typical 
velocities, around 400 m/s, are still noised at a 30% level. 
Nevertheless, one can recognize, on the wavelet-filtered 
maps, common patterns between the two time series, all 
the more easily when the scale is large. 

The next steps are now obvious: with new cameras 
with fast reading sensors (like CMOS), it is easy to in- 
crease the time sampling by a factor 10. In this case the 
k — lu filtering will be much more efficient at reducing 
the noise on the granule motion. Hence, with an increased 
time-sampling we expect to reduce the noise in two ways: 
first, by an improved efficiency of the k — uj filter and sec- 
ond by a factor \fN from the N images sharing the same 
solar signal. 

The tools developed here seem to perform quite effi- 
ciently on the numerical side and therefore allow for the 
treatment of much larger fields of view. 

Finally, we did not discuss the results in terms of solar 
turbulence. Let us mention that numbers, like the ampli- 
tude of velocity o n a 2.5 Mm scale, agr ee with previous de- 
terminations (e.g. lBrandt et allll99ll ). However, it is clear 
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that (solar) fluid mechanics should be discussed, such as 
the scale dependence of flow features. This is beyond the 
scope of this paper but will be the subject of forthcoming 
work. 



